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Abstract. We propose model of propagation of lumps of Strange Quark Matter 
(strangelets) through the atmosphere, which accounts for their apparent strong 
penetrability and normal nuclear-type sizes at the same time. The mass spec- 
trum of strangelets reaching the Earth predicted by this model agrees very well 
with the existing data on the abundance of different elemets in the Universe. 

1. Introduction 

O ' 

\& • In astrophysical literature one finds a number of phenomena which can be regarded 

as possible manifestation of the strange quark matter (SQM), extremaly interesting 
r"| ■ possibility of new stable form of matter §], §, [| ||. There were attempts to find 

lumps of SQM called strangelets in the terrestial experiments devoted to search for 
Quark Gluon Plasma but so far without success M (what can be interpreted that 
] no QGP was formed so far in these experiments j p|). The strangelets were search 

' for also in cosmic ray (CR) events [ p[ M , p"l| , [l^,]l3|, the most famous candidate 

being the so called Centauro events In CR one either deals with strangelets 
formed by some astrophysical mechanisms [[[J [l^] or one witnesses their production 
r> | proceeding in collisions of the original CR flux with atmospheric nuclei p6[ . 

■ Centauro events are good example of problem being the topic of our presenta- 

tion. They were observed very deep in the atmosphere (at ~ 500 g/cm 2 ) and if they 
would be really caused by strangelets |16| , it would require unusual penetrability 
of such objects implying their small geometrical sizes (much smaller than the typ- 
ical nuclear size). Similar conclusions were also reached in some CR experiments 
[]l2"| , Q. We present an alternative approach to the propagation of strangelets 
through the atmosphere in which we argue that they have most probably nor- 
mal nuclear-type sizes and similar propagation properties. We then demonstrate 
how such strangelets can penetrate deeply in the atmosphere producing effects as- 
cribed before to higly penetrating objects. Finally, we estimate the initial flux of 
strangelets entering atmosphere emerging from our picture and compare it to the 
known abundances of normal nuclei in Universe. 
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2. Strangelets 

Typical SQM consists of roughly equal number of u, d and s quarks and it has been 
found to be the true ground state of QCD |[[|, i.e., it is absolutely stable at high 
mass numbers A. However, any SQM produced at very early stage of the history 
of Universe would have evaporated long time ago B. On the other hand, there 
are places where SQM may still exist at present It is probably continously 

produced in neutron stars with a superdense quark surface and in quarks stars with 
a thin nucleon envelope jl5| [[J. Collisions of such objects could therefore produce 
small lumps of SQM, strangelets with 10 2 < A < 10 6 , permeating the Galaxy and 
possibly reaching also the Earth (i.e., a priori being detectable here). 

The practical measure of the stability of strangelet is provided by the so called 
separation energy dE/dA, i.e., energy, which is required to remove a single baryon 
from a strangelet. For example, if dE/dA > itin then strangelet can evaporate 
neutrons from its surface. This energy depends, among other things, on the size 
of strangelet, i.e., on its mass number A ||. There exists therefore some critical 
size A = A cr n (~ 300 — 400 depending on the parameters used) such that for 
A < A cr it strangelets are unstable against neutron emission ||. 

Anomalous massive particles, which can be interpreted as strangelets, have 
been apparently observed in three independent CR experiments: 
(i) In counter experiment devoted to study primary CR nuclei two anomalous 
events have been observed [0. They are consistent with values of charge Z ~ 14 
and of mass numbers A ~ 350 and ~ 450 and cannot be accounted for by the 
conventional background. Such values of Z and A are fully consistent with the 



(ii) The so called Price's event |l0[ with Z ~ 46 and A > 1000, regarded previously 
as possible candidate for magnetic monopole, turned out to be fully consistent with 
the above ratio Z/A for SQM @. 

(in) The so called Exotic Track with Z ~ 20 and A ~ 460 has been reported 
in [|llj. The name comes from the fact that although it was observed in emulsion 
chamber exposed to CR on balloon at the atmospheric depth of only 11.7 g/cm 2 
its arrival zenith angle of 87.4 deg means that the projectile causing that event 
traversed ~ 200 g/cm 2 of atmosphere (in contract to events (i) and (ii) where the 
corresponding depths were of the order of 5 — 15 g/cm 2 of atmosphere only). 

The Exotic Track event motivated the (balloon born emulsion chamber) JACEE 
fll2[ and Concorde aircraft |l3| experiments to search for strangelets with such long 
mean free paths. In fact, authors of |l^, [l3| suggest that the interaction mean free 
path for strangelets in atmosphere is of the order of As = 124 g/cm 2 for A = 100 
decreasing to As = 59 g/cm 2 only for A — 1000. These values are surprisingly close 
to that for protons at comparable energies (\ pr oton = 60 —70 g/cm 2 ) and are much 
bigger than that for a normal nucleus (Xnucieus — 3.8 g/cm 2 for A = 100) [^9|. It 
was then suggested in |l^, [l3| that strangelets should have geometrical radii much 
smaller than those of the ordinary nuclei or, correspondingly, much smaller inter- 
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action cross sections (in agreement with SQM interpretation of Centauro events 
mentioned before pi). 
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Fig. 1. Dependence of the rescaled radius tq [fm] (eq. ([!])): 
(a) on the ratio of s quark mass m to its chemical potential fi, 
m/fi (solid lines, read from top to bottom, correspond to fixed 
m = 100, 150, 200, 250 MeV; dashed read from top to bot- 
tom correspond to fixed /i = 100, 200, 300, 400 MeV); (6) on the 
strange quark mass m for different values of a c = 0.9, 0.6, 0.3, 0.0 
(from top to bottom, respectively). In both cases dotted line 
shows r = 1.12 fm corresponding to normal nuclear density 
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However, such expectations are not confirmed by calculations, cf. Fig. 1. 
Let us consider a strangelet visualised as Fermi gas of u, d and s quarks, with 
total mass number A, confined in spherical volume V ~ A of radius R — tq A 1 / 3 , 
where the rescaled radius rp is determined by the number density of strange matter, 



A/V — (n u + rid + n s ) /3, for which ro = (3/47rn) 



1/3 



Here n, = 



m, 



and thermodynamical potentials f2,(?7ij, Hi) are related to chemical potentials 
||. Because in our case all //, ~ 300 MeV, one can neglect the (current) masses of 
u and d quarks and keep only the s quark mass m. Taking into account the QCD 
0(a c ) corrections in fij(mj, /i,, a c ) and renormalizing them at mjv/3 = 313 MeV 
[01, one arrives at 



1/3 



r 



3tt 



2 1 



2ct c N 



/i 3 + (/i 2 — m 2 ) 



3/2 



(1) 



As one can see in Fig. la for m ~ 150 MeV and fi ~ 300 MeV (the values commonly 
accepted for SQM p], B[), the values of tq of the strangelets are comparable with 
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those for the ordinary nuclear matter. Fig. lb summarizes dependence of r$ on a c . 
As one can see QCD corrections lead to slight increase of Tq. (It is worth to say 
that at this point we differ drastically from previous works looking at strangelets 
from the QGP and Centauro production perspective [jl6| where large value of /z was 
the main reason for small values of rg and for all further claims concerning large 
penetrability of strangelets. Such large values of fi were justified there because 
strangelets were supposed to be produced from the QGP phase formed in baryon 
rich environment whereas in our case they are supposed to be stable object entering 
Earth atmosphere as a part of CR radiation.). Fig. 1 shows therefore that the 
expected decrease of the radius of strangelet is nowhere as dramatic as it has been 
estimated in Refs. [|l2], |f3|, i.e., the expected geometrical cross sections of SQM 
are not able to explain alone the occurences of anomalous events detected deeply 
in the atmosphere. 

3. Possible scenario of propagation of strangelets through the 
atmosphere 

It does not mean, however, that some form of SQM does not reach Earth and pen- 
etrates deep in the atmosphere to be finally registered. The apparent contradiction 
between its "normal" size and strong penetrability can be resolved in very simple 
manner. In what follows we shall propose the respective (speculative but plausible) 
scenario, which can be summarized in the following way: Strangelets reaching so 
deeply into atmosphere are formed in many successive interactions with air nuclei of 
much heavier lumps of SQM entering our atmosphere. In this way one can accomo- 
date both the most probably "normal" mean free paths for successive interactions 
and final large penetrating depth. Such scenario is fully consistent with all present 
and proposed experiments and could be additionally checked only by measuring 
products of the intermediate collisions, which so far is impossible. 

To provide numerical estimations we shall limit ourselves to two extreme pic- 
tures of collisions of strangelet of mass number A with air nucleus target of mass 
number A t : 

1. "Standard model": all quarks of A t which are located in the geometrical inter- 
section of two colliding projectiles are involved and it is assumed that each quark 
from the target interacts with only one quark from the strangelet; i.e., during inter- 
action 3A t quarks from strangelet can be used up and after it the mass number of 
strangelet is diminished to the value of A — A t a.t most. In this case simple analytical 
estimation is possible leading to the total penetration depth of strangelet 

*-i^(^r(*-^)'(-^)-^(ir- <*> 

(where Ajv,4 t denotes the mean free path for NA t collisions). 

2. "Tube-like model": All quarks from both nuclei which are in their geometrical 
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intersection region participate in the collision. This is an analogue of the so called 
tube model (TM) used sometimes in nuclear collisions and after collision the atomic 
mass of strangelet diminishes to the value equal to A — A 1 / 3 • A^ 3 (this variant 
should be regarded as a limiting possibility providing estimation of the maximal 
possible destruction of the quarks in the strangelet) . 




Fig. 2. (a) An example of the predicted decrease of actual 
size of strangelet A with depth h of atmosphere traversed for two 
different initial sizes: Aq = 1000 and 2000 and for two scenar- 
ios of interaction: standard (solid and dotted lines) and tube-like 
(dashed and dotted lines). In both scenarios dotted lines corre- 
spond to A < A cr i t . (b) Atmospheric length A after which initial 
strangelet reaches its critical dimension, A = A crit (= 320) drawn 
as a function of its initial mass number A for standard (solid 
line) and tube-like (dashed line) interaction. Consecutive squares 
indicate points where (for A > 600) A /A crit = k, (k = 2, 3, . . .). 

In the numerical calculations we take A t = 14.5 for the atmosphere (assuming 
25% of oxygen and 75% of nitrogen) and neglect the influence of the interstellar gas. 
As a result one obtains the mass number of strangelet registered at depth h, A(h), as 
a function of h as shown in Fig. 2a. As one can see, bigger initial strangelets (i.e., 
with higher mass number Aq) can penetrate much more deeply into atmosphere 
untill A(h) exceeds critical A crit , after which point they just evaporate by the 
emission of neutrons. Notice much stronger degradation of the initial strangelets 
for the case of the tube-like model interactions with air nuclei. As we have already 
mentioned, our scheme displays two most extreme situations in order to provide 
limits for all further possible estimations. Because, for simplicity, we tacidly assume 
that air nucleus destroys totally only the corresponding (equal to it) part of the 



Propagation of Strange Quark Matter. 



6 



incoming strangelet, our estimation provides therefore only a lower limit of what 
should be expected in more detailed calculations. More refined estimations should 
take into account also the surface tension effects and the possible dependence of r 
and cross section on the actual mass number of strangelet A(h), which translates 
immediately into dependence on h. However, their effect would be of secondary 
importance only and at present experimental situation such detailed and tedious 
estimations would be simply premature. 




Fig. 3. (a) The estimation of the expected flux of strangelets on 
the border of atmosphere, N(Aq), as a function of their mass num- 
ber as obtained from standard (full symbols; solid line ~ Aq 7 ' 5 
and dotted line ~ exp (— Aq/130) ) and tube-like (empty symbols; 
dashed line ~ A^ 2 ' s ) models. The respective points are taken 
from: Saito §L E T @, Centauro @, JACEE Concorde 
jl3| and Price |lf|. See text for further details, (b) Comparison of 
N(Aq) from (a) with the abundance of elements in the Universe 
after pifl . 

Suppose now that the energy per baryon in strange matter is e = 919 MeV 
and that number densities corresponding to nuclear matter are n = (ffO MeV) 3 
||. In such situation, for A < 1 fOO, E/A exceeds already tojv but strangelet does 
not emit neutrons yet and starts to do so only for A < 320, at which point dE/dA 
exceeds || . Below this limit strangelet decays rapidly by evaporating neutrons. 
In view of these remarks it is remarkable that all possible candidates for SQM have 
mass numbers near or slightly exceeding A cr u, namely: A — 350 and 450 in 
A = 460 in |lf| and A = 1000 in jl8| and it is argued that Centauro event contains 
probably ~ 200 baryons [|f4| . 

Fig. 2b shows atmospheric length traversed after which strangelet mass num- 
ber A becomes critical, A — A crit (starting from different mass numbers A of the 
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initial strangelets) . From this figure one can read off that strangelets which are 
observed at depth 200 g/cm 2 should originate from strangelets of mass number 
Aq = 900 at the top of the atmosphere whereas Centauro events observed at the 
mountain altitudes would require original strangelet of Ao = 1800. 

Let us now estimate flux of strangelets reaching our atmosphere starting with 
the experimental data taken at different atmospheric depths, cf. Fig. 3a. The 
experimental data Saito, ET and Centauro on measured fluxes on different atmo- 
spheric depths (which can be interpreted in terms od SQM) are taken from 
p8[ and |l4f| , respectively. The corresponding upper limits (no strangelets found 
so far), JACEE and Concorde, are from H] and Jf3|. Notice that Price data |§ 
(assuming that what was observed was indeed a strangelet) favour standard model. 
In terms of fits (for 3 points: Saito, ET, Centauro) one gets: ~ Aq 7 ' 5 for standard 
model (full line, dotted line corresponds to ~ exp (— Aq/130) in this model) and 
~ Aq 2 ' 8 for tube- like model (dashed line). The choice of the power- like or exponen- 
tial form for standard model was dictated by the analogy to nuclear fragmentation 
and the expectation that decay (fragmentation) of a strange star after its collision 
will result in the production of strangelets with similar distribution of mass (2(j) . It 
is interesting to note that essentially the same exponent, A^ 7 ' 5 , as obtained in the 
standard model fit is observed also for occurence of normal nuclei in the Universe 
@, cf. Fig. 3b. 



4. Summary and conclusions 

We summarize by stating that most probable the geometrical cross sections of 
strangelets are not dramatically different from those for the ordinary nuclear matter 
and cannot therefore explain their apparent very high penetrability through the 
atmosphere. Instead we propose to interprete such a penetrability of strangelets 
(already discovered or to be yet observed) as indication of the existence of very 
heavy lumps of SQM entering our atmosphere, which are then decreased in size 
during their consecutive collisions with air nuclei (i.e., their original mass number 
Aq is reduced until A — A cr n) and finally decay be the evaporation of neutrons. 
Assuming that mass distribution of strangelets comprising this initial flux follows 
power low Aq we estimate that 7 > 2 is consistent with present experimental data. 
Our result is consistent with other estimations of the flux of strangelets at different 
atmospheric depths. In particular it is interesting to note that the best fit to data 
representing this dependence comes from our standard model and it agrees very 
well with the existing data on chemical composition of normal nuclei supporting 
therefore our hypothesis that geometrical sizes of strangelets are not much different 
from those of normal nuclei of the same A. 
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